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Abstract. The contributions of three different types of driving terms are included in the estimation of 
the pd — > pdr) reaction at low energies. Near threshold, it is predicted that a two-step model involving 
an intermediate pion should be the most important but, as the energy approaches the threshold for rj 
production in the free nucleon-nucleon reaction, a pick-up mechanism with a spectator proton should 
become dominant. The total cross sections are underestimated by about a factor of two compared to 
experimental data but the discrepancies in the angular distributions are more serious, with no sign in the 
data for the peaks corresponding to the pick-up diagram. 

PACS. 25.10.+S Nuclear reactions involving few-nucleon systems - 13.60.Le Meson production - 14.40.Aq 
pi, K, and eta mesons 



1 Introduction 

The great interest in the production of heavy mesons near 
threshold started with measurements of the pd — > 3 He?7 
reaction, which showed a surprisingly strong cross section 
with an anomalous energy dependence pQ . An impulse ap- 
proximation description of the process, where the produc- 
tion takes place in nucleon-nucleon scattering, with the 
other nucleon in the deuteron being essentially a spec- 
tator, greatly underestimates the observed rate because 
of the high momentum components required in the nu- 
clear wave functions 2j. To share the large momentum 
transfer between the nucleons, a mechanism was proposed 
whereby a pion was produced in an initial pp — > dir + re- 
action to be followed by the production of the observed 
meson through a secondary ir + n — > i]p process ,3 . Near 
threshold, the kinematics are favourable for the final pro- 
ton and deuteron to stick to form the observed 3 He. A 
quantum-mechanical evaluation of this suggestion 0] led 
to a cross section that was only about a factor of two 
lower than a precise measurement of the reaction rate |S] . 
To understand the situation further, it would be helpful 
to look at other related final states. 

The total cross section for the pd — > pdr\ reaction, 
where a final 3 He is not formed, was measured at two en- 
ergies very close to threshold at Saclay [H] . A preliminary 
evaluation of the two-step model with an intermediate 
pion, that was successful in the description of the 3 He?7 
final state, gave very promising results when compared to 
these data 0. Differential as well as total cross sections 
for this reaction at higher energies have recently become 
available from Uppsala 8 , to complement those obtained 



by the same group for the pd — + 3 He r\ reaction |jjS]. The 
time therefore seems opportune to make a further theo- 
retical investigation of the pdr] final state. 

The kinematics of the Saclay experiment j^j are very 
close to those where a 3 He emerges and so it is not sur- 
prising that in this region the two-step model describes 
the process well, as it did for the 3 Hery final state 0]. On 
the other hand, as one approaches the threshold for r\ pro- 
duction in free nucleon-nucleon collisions, one expects a 
pick-up diagram, corresponding to a quasi-free pn — > dr\ 
production on the neutron in the target, to dominate. Be- 
tween these two extremes there is also the possibility of 
a contribution from the impulse approximation diagram 
with a quasi-free pN — > pN r\ on a proton or neutron 
that is bound both initially and finally in the deuteron. 
Away from threshold, these terms may not be universally 
suppressed, as they are for 3 Hery, because of the greater 
flexibility in the pdrj kinematics; in certain parts of phase 
space, the momentum transfers between initial and final 
deuterons are minimised. 

The three types of contributions are described further 
in sect. 2, where the general kinematics are discussed. The 
following three sections are devoted to the evaluation of 
the cross section distributions for the individual terms. 
Uncertainties in the phases of the amplitudes leads us 
to neglecting interferences. This may not be too danger- 
ous because the various models tend to populate different 
parts of what is a three-body phase space. More worrying 
is that we neglect also all final-state interactions fsi. Now, 
unlike the r\ 3 He case [5], there is no simple prescription for 
the inclusion of such effects when there is an fsi between 
more than one pair of particles. Experimentally a thresh- 
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old enhancement is seen only in the r\ d invariant mass 
distribution, whereas that for pd looks like phase space, 
despite the presence of the strong interaction which could 
lead to the formation of the 3 He [H]- The results for the 
total cross sections, and both angular and invariant mass 
distributions, are compared with experiment in sect. 6. 
The overall production rate is underestimated by about 
a factor of two, which is rather similar to that found for 
the pd — ► 3 He?7 reaction within a similar theoretical ap- 
proach. However, the angular distributions of both the 
proton and deuteron show discrepancies with respect to 
the data, with the latter showing no signs of the forward 
deuteron and backward proton peaks corresponding to the 
spectator proton of the pick-up mechanism. Our conclu- 
sions are drawn in sect. 7. 



2 Defining the problem 

The three classes of diagram that we evaluate for the 
pd — ► pdr] reaction are illustrated in Fig. ^ It is gen- 
erally assumed that the neutron-exchange diagram (a) 
dominates the reaction above the free NN threshold and, 
indeed, it is this hypothesis that is the basis of the extrac- 
tion of the quasi-free pn — > drj cross section on a moving 
neutron from pd — > pdr] data ^Oj ■ Below threshold, the 
presence of a spectator proton in this case will bias the 
distributions to low d r] invariant masses. 

The impulse approximation diagram of Fig. \\]p is of- 
ten used to model coherent reactions on the deuteron at 
high energies, where the momentum transfer to the final 
deuteron can be small. There is, of course, a second term 
where the r] production takes place on the neutron and 
some theoretical model is required to deduce the rela- 
tive phases of the spin-isospin input amplitudes needed 
to make reliable estimates here [111 . A final-state interac- 
tion between the proton and deuteron to form an 3 He for 
either the diagrams in Fig. fjji,b would lead to the triangle 
graph, whose contribution has been shown to be small for 
pd — > 3 He 77 FJ]- That they are not necessarily negligible 
here is due to the fact that the final dp pair does not have 
to emerge with low excitation energy. 

The differential cross section for the pd — > pdr] re- 
action is determined by the matrix element M. through 




Fig. 1. Three classes of diagrams relevant for the pd — ► pdr] 
reaction, (a) The pick-up term which dominates above the 
NN threshold, (b) The impulse (triangle) diagram. There is 
a similar contribution with the internal proton and neutron 
interchanged, (c) The two-step model with an intermediate 
pion. The contribution where the 7r + is replaced by a ir° is 
related to this through isospin invariance. 



(2.1) 

where the sum is over the spin projections in the initial 
and final states. Here Pf — Pj = s is the square of the 
total cm energy and {Ej,p^) are the energy and three- 
momentum of the reaction products. The incident flux 
factor in the centre of mass involves the initial proton cm 
momentum p. 

Experimental distributions have been presented in terms 
of the angles of the final particles and the invariant masses 
of pairs of such particles [Hj ■ For this purpose it is conve- 



nient to reduce the phase-space factors to: 




p d dS7 d p*dQ* dm. 



(2.2) 

where the final deuteron energy and momentum (Ed,Pd), 
as well as the angles (Id, are evaluated in the overall cm 
system while those of the r] are evaluated in the rjp rest 
frame, where the invariant mass is given by 

ml p = s + m d -2</~sE d . (2.3) 
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Analogous formulae follow immediately for the other two 
combinations of variables. 



3 The pick-up contribution 

Though there is an enhancement at threshold, the energy 
dependence of the pn — > d i] total cross section is broadly 
consistent with s-wave production up to an excess energy 
of at least 60 MeV • At threshold there is only one cm 
amplitude, which can be written as 



M{pn — > dn) = ul c [Gy/s^p n ■ e d ] u p , 



(3.1) 



where p n is the initial momentum and the cm energy ^/s pn 
arises from the reduction from a relativistic form. Here u p 
and u„ e are proton and charge-conjugate neutron Pauli 
spinors respectively. 

In terms of the amplitude G, the total pn — > dn cross 
section is 



<jt {p n —> drj) 



1 

8^- 



(3.2) 



where the p* are evaluated in the pn cm frame. The ex- 
perimental data are then consistent with a value of \G\ 2 = 
(0.046 ± 0.011) fm 4 [Ej. 

The spin-averaged amplitude squared for the pick-up 
term of Fig. reduces to: 

\\M(pd P d v )\ 2 

= \[{2irf2m d ](AE pPp ) 2 \G\ 2 {f> s (q) 2 + ^n(q) 2 } (3.3) 

Here the deuteron wave functions <fis(<j) and <Pd(<1) are 
normalised by 



&s(q) 2 + ^D(q) 2 } = 1. 



They are evaluated at a momentum-squared of 

2 

2 ^^ri /I \ 2 2 

1 =^t{-W-PvV +P P ±> 



(3.4) 



(3.5) 



where the Lorentz boost has been approximated in the low 
Fermi momentum limit and {p p \\ ,P P ±) are the components 
of the final proton momentum p p parallel and perpendic- 
ular to the momentum p of the incident deuteron, which 
has energy E dl . 



4 The triangle diagram 

The triangle diagram of Fig.^i requires as input the am- 
plitudes for the sub-reactions NN — > NNi]. The thresh- 
old cm amplitudes for 1=1 and 1 = are respectively: 



Mi(NN -> NNrj) 



Wi, s rjj p ■ e 
M (NN NNr,) = \ W 0< t V ■ e/ V 



t 



Xf-Xi, (4- 1 ) 
(4.2) 



where \ and <f> represent the isospin-1 and isospin-0 con- 
figurations of the NN states, with e and r\ corresponding 
to the spin-1 and spin-0 combinations. 

The spin-averaged total r\ production cross sections 

are 



<t{pp^ppv) = 



m N 



512ir 2 p NN s NN (i + 2m N / mri ) 1/2 



*Q 2 nn v \Wi,s\ 2 , 



o~(pn—>pnn) 



m N 



1024:Tr 2 p NN s NN (i + 2m N /m v ) 1/2 
Q 2 nn v x [|W M | 2 + |yM 2 ] , (4.3) 

where pnn and sjvjv are the momentum and square of the 
total energy in the cm system, mjy the nucleon mass, and 
Qnnti the excess energy in the final system |llj . 

In terms of an initial deuteron momentum p, a fi- 
nal deuteron momentum p d , and a momentum transfer 
Q = p d + p, the triangle graph is evaluated in the stan- 
dard way by putting the spectator nucleon on-shell at an 
average energy Em{\p)- The spin-averaged matrix ele- 
ment becomes 



1 ~Y,\M(pd^pdr 1 ) \ 2 =±- 



spin 



m d 



64 \AE N (p/2) 



[|3Wi, a + W , t \ 2 + S 2 Q (\Q) + |5i(|Q)] 

+\3Wi, s -W ,t\ 2 IS 2 M (±Q)] . (4.4) 

The spherical, quadrupole, magnetic and axial form 
factors appearing here are defined by 

S S (Q) = S a (Q) + S b (Q) , 

S Q (Q) = 2S c (Q)-S d (Q)/V2 

S M (Q) = S a (Q) - \S b (Q) + S C (Q)/V2 + \S d {Q) , 

Sa(Q) = S a (Q) - |5 6 (Q) - V2S C (Q) - S d (Q) , (4.5) 

where the terms on the right hand side represent integrals 
over the reduced deuteron wave functions u(r) and w(r): 



S a (Q) = 


/ Jo (Or) 
Jo 


[it(r)] 2 dr , 




S b (Q) = 


poo 

/ Jo (Or) 
Jo 


[w(r)] 2 dr , 




Sc(Q) = 


/•OO 

/ MQr) 
Jo 


u(r)w(r) dr , 




S d (Q) = 


/>oc 

/ J 2 (0r) 
Jo 


[w(r)] 2 dr . 


(4.6) 



Experimental data on n production in pp |12| and pn |13j 
collisions above the fsi regions show that | VKi . s | 2 m 1.9 x 
10 8 fib and |Wo,t | 2 « 23 x 10 8 /ib, though the relative 
phase of these two amplitudes is not determined. In a 
meson-exchange approach, it is suggested that p-cxchange 
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is the largest term |1 lj . If further we take the two ampli- 
tudes to be in phase then we find that 



|3W M 
|3W M 



W ,t| 2 
W ,t| 2 



80 x I0 8 pib , 
0.4 x 10 8 /ib . 



(4.7) 



5 The two-step model 



The contribution of the two-step model with an interme- 
diate pion can be estimated using the same techniques 
as for the more complicated pd — > 3 He r\ reaction, where 
there arc two integration loops over the unobserved Fermi 
momenta 0] . In terms of the internal momenta marked in 
Fig. lc, the pd — > pdr] matrix element becomes 



M 



(2tt) 3 / 2 



d 3 q' M(ir + n — > r/p) x 

int 

M(pp^ 7r+d)#(g') , (5.1) 



where the sum runs over the internal spin indices. 

The pion propagator between the matrix elements for 
the production and conversion of the pion in eq. 15.111 has 
been approximated by its positive energy pole. The dif- 
ference between the external and internal energies, AE — 
V^s — E int , depends upon the Fermi momentum q and, 
following ref. 0], we retain only linear terms in this inte- 
gration variable: 



AE = yfs 



s - E d { Pd ) - 
AEo + V-q, 



E n {- 



q) -E n {\p-q 



where the mean energy defect 

AE Q = - E d ( Pd ) - E n (-ip) - E v (\p - p d ) 



-Pd) 
(5.2) 



(5.3) 



and the relativistic relative velocity between the pion and 
neutron 



V 



E n + Ett 
2E n E-rr 



1 



(5.4) 



q occurs is in the denominator of the propagator. By defin- 
ing a Lorentz transformed velocity with components 



V' ± = V± and V'» = 7V11 



(5.5) 



the whole integrand can be written in terms of q' . 

After decomposing the deuteron wave function into its 
I ~ 0, 2 components, we are left with integrals of the form 



Ml')Ti(q') 



AE + V' ■ q' + ie 



dt e t ( AE o+' ie ) t 
q' 2 dq'<p e (q') I AQ q , e uV ' ^ T e (q>) 



S t {AE ,\V'\)Ti(V) 



(5.6) 



where the S- and Z?-state form factors involve integrals 
over the configuration-space deuteron wave functions 



Si(AE ,\V'\) = 



(27T) 



3/2 



V 



dtt 



Mt) (5.7) 



with uj = AE /\V'\. 

The s-wave nature of the 5n(1535) resonance that 



dominates low energy 77 production in 



VP, 



leads 



to largely isotropic production at low energies. We can 
then take the matrix element to be proportional to the 
corresponding spin-non-flip amplitude 



A4(ir + n — > r/p) = 



47T, 



VP 



VP) 



(5.8) 



For simplicity, of the six invariant pp — > ir + d ampli- 
tudes we shall keep only the largest in our energy region, 
for which 



M(pp 



-d) 



-V2A(e d -p n ), 



'pp 



(5.9) 



where 4> pp and e d represent the spin-zero and spin-one 
initial and final TV TV states. This is then related to the cm 
cross section through 



|^| 2 =8(27r) 



m 2 p w df2 



ir+d) 



(5.10) 



do not depend on the proton or 77 momenta in the final 
state. 

Both the argument of the deuteron momentum space 
wave function ^{q') and the integration variable are the 
Fermi momentum q' , Lorentz contracted in the beam di- 
rection; q' ± = q ± , qj| = q^/j, with 7 = E d (p)/m d w 

E n {~\p)jm n . 

The linearisation is valid only if relatively small Fermi 
momenta are required to allow the two steps to proceed al- 
most on-shell, which means that the energy defect should 
be small. At the threshold for n production AEo m — 15 MeV 
which illustrates the kinematic miracle first noted for the 
pd — > 3 He?7 reaction 0. We then neglect the dependence 
of the individual pp — > ir + d and tt + u — * i]p amplitudes 
upon the Fermi momentum so that the only place where 



In addition to the tt + diagram of fig. lc, there is also 
the possibility of tt° propagation between the two steps. 
This can be taken into account by simply multiplying the 
matrix element M of eq. (5.1 ) by an isospin factor of 3/2. 
Using eq. I|2.2|l . we then arrive at an expression for the 
distribution of the pd — > pdrj cross section in terms of 
the final deuteron angles and rjp invariant mass: 



d a 



9 



Sr/p TUdPdPr, 



df2 d dm np 4(2tt) 4 spEl 



(\SsY 



\St 



pp 

m 2 Pit 



p p da 



an 



(PP ~ ¥ TT + d) 



Pr, 



r)p) 



(5.11) 



where we have taken advantage of the isotropy of the 
n + n — » np differential cross section to integrate over the 
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77 angles and to write the result in terms of the total cross 
section. This option is not open for the other differential 
distributions, where extra non-trivial integrations have to 
be performed. 

An effective range description of the low-energy n + n — > 
rjp total cross section is 



Pn 



(7r + n 



VP) 



2.76|a| 2 



1 — iap* + 



^ r oap^ 2 \ 



(5.12) 



The data base has not improved significantly since we took 
the values a = (0.476 + 0.279i)fm and r = (-3.16 - 
0.13i)fm gj. 

The pp — » 7r + d cross sections are obtainable from the 
SAID analysis ^3] but, because the intermediate pion in 
Fig.^L is non-physical, some prescription is needed to ex- 
trapolate from the experimental data. We have assumed in 
the applications that the amplitudes are unchanged when 
the cm production angle is kept fixed. 



6 Comparison with experiment 



c 













A'' 






■7 / 












1 


50 

Q (MeV) 


100 



Fig. 2. Total cross section for the reaction pd — > pdr/. The 
experimental data are taken from Saclay (solid circles) and 
CELSIUS (open circles) The dotted curve represents the 
predictions of the triangle diagram (impulse approximation), 
the chain curve the pick-up contribution, and the dashed curve 
that of the two-step model. The incoherent sum of these three 
terms is shown by the solid curve. 

The predictions of the three driving terms of Fig. ^ are 
compared in Fig. [21 with the existing experimental data |SJ 
[S] . It is immediately apparent that the impulse approxima- 
tion of the triangle graph in Fig. lies about two orders 
of magnitude below the experimental data. The ambiguity 
in the relative phase of the amplitudes Wi s and Wo,* is 



therefore largely irrelevant and the impulse approximation 
term will now be dropped from further discussion. 

In the near-threshold region of low Q, the two-step 
contribution dominates that of the pick-up by an order 
of magnitude. The kinematics here are somewhat similar 
to those of low-energy pd -> 3 He 77 reaction where the 
pick-up term with a proton-deuteron final-state interac- 
tion similarly underpredicts the experimental data |15ll2] . 
However, the two-step term levels off at Q « 80 MeV and 
eventually decreases at higher energies due to a combina- 
tion of factors. The region where the intermediate pion is 
close to being physical becomes a smaller fraction of the 
allowed phase space but also the pp — > dir + cross section 
falls steeply above the A region. Since the pick-up term 
must approach that for quasi-free pn — > drj above the NN 
threshold, the two contributions must cross before then. 
In our estimates, this occurs at Q « 95 MeV. 

The incoherent sum of the three contributions is also 
shown in Fig. [21 Though the shape is very similar to that 
of the experimental data, these lie about a factor of two 
above the curve at high energies and a little bit more near 
threshold. This parallels a similar underprediction of the 
pd — > 3 He 77 reaction in a two-step model that includes 
the effects of the 77 3 He fsi 0] . The fractions of phase space 
where the pd or 77<i final-state interactions might enhance 
the pdr] channel are relatively small at the higher energies. 

More information can be derived from the differen- 
tial distributions that were measured at T p — 1032 MeV 
(Q = 72.3 MeV) jHj- The only significant enhancement in 
the three invariant mass distributions of Fig.|3|is that near 
the rjd threshold, which is consistent with the large scat- 
tering length suggested by the low energy pn — ► ^77 total 
cross section data ^UJ- I n contrast, despite the existence 
of the nearby 3 He bound-state pole, there is no sign of 
any pd fsi. The two-step model leads to only minor dis- 
tortions of the invariant mass phase spaces, pushing m m 
to slightly larger values, whereas the pick-up contribution 
preferentially populates low masses which is reflected 
kinematically as a slight increase at higher values of m p d- 

Apart from the overall strength being too low by a fac- 
tor of two, the sums of the two contributions give plausi- 
ble descriptions of the invariant mass distributions but the 
same cannot be said for the angular distributions shown 
in Fig. 0] In the pick-up contribution of Fig. QJi, the fi- 
nal proton is a spectator and the transverse momentum 
is governed by the Fermi momentum components in the 
deuteron. This automatically gives a sharp peak for a cm 
proton angle close to the backward direction and this has 
a slightly weaker kinematic reflection around the forward 
deuteron direction. Given that the estimate within the 
pick-up model has relatively few uncertainties, the dis- 
crepancy with the experimental data of Ref. [5| is partic- 
ularly significant. 

There are far more uncertainties in the evaluation of 
the two-step model since the intermediate pion is gener- 
ally off its mass shell and this effects the kinematics of 
the pp — > dir + reaction. However, the deuteron angular 
distribution of Fig. does suggest that the models need 
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Fig. 3. Distributions in the invariant masses of thepd — > pdr; 
reaction at Q = 72.3 MeV. The panels (a), (b), and (c) refer 
respectively to the prj, pd, and drj cases. The experimental data 
from Ref. [Sj are compared to the predictions of the pick-up 
term (chain curve), the two-step model (dashed curve), and 
their incoherent sum (solid curve). 



o 

COST? 



Fig. 4. Centre-of-mass angular distributions with respect to 
the incident proton direction of the pd — > pdr) reaction at 
Q — 72.3 MeV. The panels (a), (b), and (c) refer respectively 
to outgoing d, rj, and p. The experimental data from Ref. [B] are 
compared to the predictions of the pick-up term (chain curve) , 
the two-step model (dashed curve), and their incoherent sum 
(solid curve). 



to be supplemented at large angles and this might be the 
reason for the underestimation of the total cross section. 



7 Conclusions 

We have estimated the contributions of three different 
models to the total and partial cross sections of the pd — > 
pdrj reaction below the 77-production threshold in nucleon- 
nuclcon collisions. The impulse approximation turns out 
to be largely negligible compared to the other two terms 
of which the two-step model is shown to be dominant in 
the near-threshold region. However, although some of our 
kinematic approximations may start to break down before 



one reaches the iieepn threshold (Q w 192 MeV), the esti- 
mates suggest that the two-step mechanism would provide 
only a very small correction to the pick-up interpretation 
of the CELSIUS quasi-free pn — + dr\ data |TU]. 

We have neglected the final-state interactions which 
should distort the lower edges of the invariant mass dis- 
tributions shown in Fig. [3] while having a smaller effect 
on the total cross section away from the threshold region. 
In fact the only fsi clearly seen in the experimental data 
is that associated with the rjd channel [Sj. Now the to- 
tal cross section at the higher energies is underestimated 
by a factor of two, which is a very similar factor to that 
found for the near-threshold two-body pd — > 3 Hery re- 
action interpreted in the same two-step approach, after 
the inclusion of the necessary fsi effects 0]. It is, how- 
ever, intriguing to note that, although there is no sign of 
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the strong pd fsi in the experimental data of Fig. Eb, the 
normalisation of the total pd — > pdr\ total cross section 
is predicted successfully from the pd —> 3 He r\ data using 
the final-state extrapolation theorem ^H] without imple- 
menting a full dynamical model. 

Apart from the fsi regions, it is hard to draw firm 
conclusions on the models from the invariant mass dis- 
tributions of Fig. [3J There is far more information to 
be gleaned from the angular distributions of Fig. 0] The 
models underpredict the data in the backward deuteron 
hemisphere and, more critically, there is no sign in the 
data for the sharp peaks for forward-going deuterons and 
backward-going protons being produced in the pick-up 
process. Though these might be softened somewhat by 
multiple scatterings or fsi, it should also be noted that if 
cither the deuteron or the proton is lost down the beam 
pipe then the event is not registered [H]. 

If indeed the two-step model is the dominant mecha- 
nism well below the nucleon-nucleon threshold, it would 
have consequences for the pd — > ppnrj four-body final 
state, where one could expect to see a strong pn fsi. Data 
on this reaction were taken at CELSIUS simultaneously 
with those on the other channels, but their detailed anal- 
ysis is not yet complete f^Z]. Data also exist from COSY 
on the analogous sub-threshold pd^ K + A d reaction JHj 
and it would seem likely that the corresponding two-step 
model should play an important role there as well. 

We are much indebted to our colleagues in Uppsala, especially 
J. Zlomahczuk, for discussions regarding the ^-production pro- 
gramme at CELSIUS. One of the authors (CW) is grateful for 
financial support and hospitality from both the Department of 
Radiation Sciences and the The Svedberg Laboratory of the 
University of Uppsala. Comments by Y. Uzikov regarding the 
impulse approximation term proved helpful. 
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